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The Lagrangian description of fluid flow in relativistic cosmology is ex- 

o 

tended to the case of flow accelerated by pressure. In the description, the 
entropy and the vorticity are obtained exactly for the barotropic equation of 
. state. In order to determine the metric, the Einstein equation is solved per- 

turbatively, when metric fluctuations are small but entropy inhomogeneities 
are large. Thus, the present formalism is applicable to the case when the 
, inhomogeneities are small in the large scale but locally nonlinear. 
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I. INTRODUCTION 



o 

H . 

Understanding the evolution of fluids in the expanding universe is of great interest in 
the cosmology. Here, we develop the relativistic version of the Lagrangian description of 
fluid flow for the following reason: First, let us make a comparison between the Eulerian 



description and the Lagrangian one. In the Eulerian approach, all variables are expanded 
in series. Consequently, the density contrast should be small. On the other hand, the 
Lagrangian framework is based on the description along the fluid flow, so that we can solve 
exactly the continuity equation for the density. This is a great advantage, since we can tackle 
with the dynamics at the nonlinear stage. Next, we compare the Newtonian treatment with 
the general relativistic one. The Newtonian treatment is often used as a good approximation 
for the region l/L <C 1, where / is the length scale of fluctuations of fluids and L corresponds 



*Electronic address: asada@phys.hirosaki-u.ac.jp 



1 



to the Hubble radius. Thus, the treatment is restricted within small scales, though it enables 
us to understand its results quite intuitively. On the other hand, there are no restrictions on 
scales in the general relativistic treatment. Hence, the relativistic version of the Lagrangian 
description is most useful for studying highly nonlinear region in the expanding universe 
up to the caustic formation. For dust, the Lagrangian description is formulated generally 
|l6|,|2]fj| . Now, we extend it to the case of fluid with pressure. This case include the radiation 
dominated era in the early universe and collapsing region in which the velocity dispersion 
grows significantly at the late stage. 

A key idea in the Lagrangian description of the universe is as follows We illustrate 
it in the Newtonian cosmology for simplicity. The density is exactly obtained along the fluid 
flow. The Poisson equation for the cosmological Newtonian potential <fi is 



where 6 denotes the density contrast. In the cosmological situation, can be safely consid- 
ered as small, even if the density contrast blows up in the small scale. Actually, this occurs 
when l/L goes to zero, namely in the small scale. We can expect that the idea works well 
also in the relativistic case, if <fi is replaced by the metric. The fluid flow approach has been 
discussed by several authors However, their treatment is not satisfactory, since their 

formulations are based on the fluid flow but they solve their equations perturbatively by 
splitting the flow into the background and perturbed parts. Hence, their approach is still 
restricted within the small density contrast. 

In section 2, we consider the perfect fluid. It is shown that the entropy and the vorticity 
are determined exactly based on the Lagrange condition. Under this condition, section 3 
presents a perturbative Lagrangian approach. Conclusion is given in section 4. Greek indices 
run from to 3, and Latin indices from 1 to 3. We use the unit, c = 1. 



A0 = 4nG(p- p b ), 



(1.1) 



where pb denotes a density in a background universe. This is estimated as 




(1.2) 
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II. LAGRANGIAN DESCRIPTION 

Let us consider a universe filled with the perfect fluid, whose energy momentum tensor 
is written as 

T Mi/ = (e + P)uV + Pg^. (2.1) 
The conservation law becomes 

(£u") ;M + iVV = 0, (2.2) 
(e + PKX + ^Vtf = 0, (2.3) 

where we defined the projection tensor as 

Y u =K + ^u v . (2.4) 

Here, we assume the barotropic equation of state as P = P(e), so that we can introduce 
the entropy and the enthalpy respectively as [HI, IT| 



r de 

s = expj j—p, (2-5) 
r dP 

Then, the conservation law is re-expressed as 

= 0, (2.7) 

u^u v + {hxh) >vl ;^0. (2.8) 

Furthermore, we define the vorticity as 

^ = \ e ^ UaU ^ (2.9) 

where denotes the complete anti-symmetric tensor with e 0123 = 1/^—g and g = 

det((7 AiJy ). Then, the Beltrami's equation for the vorticity is written as 

« %. = (%) **. (2.10) 

v s ' > v h v h ' > s 



To this point the treatment is fully covariant. In the following, we introduce the La- 
grangian coordinate, 



z" = (t,s*), 



(2.11) 



where r is the proper time and x % is constant along the fluid flow. This gives us the 
Lagrangian condition (e.g. |12|), in which the matter 4- velocity take components of 



vJ 1 = (1,0,0,0). 



(2.12) 



Under this condition, we have g 00 = —1 and w M = (—l,g 0i ). In the Lagrangian description, 
the entropy conservation equation Q2.7|) is simply 



[s^) fi = 0. 



(2.13) 



Therefore, 



s(x, t) 



-s(x, t ). 



The relativistic Beltrami equation ( |2.1Up also becomes simply 



(2.14) 



(2.15) 



which is integrated to give 



hu l hu) 1 



(2.16) 



to 



This is also expressed as 



huj l (x, t 



o 



The u° component is not independent of u l : From uo^u^ = 0, we obtain 



(2.17) 



u° = g 0i uj l . 



(2.18) 



The result Eq. ( |2.16p tells us that the vorticity is coupled to the entropy enhancement 
and vice versa. In particular, if the vorticity does not vanish exactly at an initial time, the 



vorticity will blow up as the entropy grows larger and larger (i.e. in the collapsing region), 
even if it has only the decaying mode in the linear perturbation theory. It should also be 
emphasized that our results Eqs. ( [2.14| ) and (|2.16| ) in the fully general relativistic treatment 



precisely correspond to those in the Newtonian case. 

III. PERTURB ATIVE APPROACH 

In the previous section, we solved exactly the equations for the entropy and the vorticity. 
The results Eqs. ( p,14|) and (|2.17|) show that s and u 1 are completely written in terms of the 



determinant of the metric tensor and their initial values. Here, we shall obtain perturbatively 
the metric at the linear order. 

The Einstein equation is decomposed with respect to the fluid flow: 

G^u u = 8nGe, (3.1) 

G F uY Q = 0, (3.2) 
G^-fa-fp = 8nGP la(3 . (3.3) 

The Euler equation fl2.8|) is rewritten as 

(/iO,o + ^ = 0. (3.4) 

We assume that the background is spatially flat Friedmann-Lemaitre-Robertson- Walker 
(FLRW) universe. The extension to the spatially non-flat case must be a straightforward 
task. The perturbed metric is decomposed into 

got = B ti (x) +h(x), (3.5) 



9ij = a 2 (5ij + 2H L 8 i:j + 2H Tjij + (h id + h j4 ) + 2H i:j ). 



where B, Hi, and Ht are scalar mode quantities, hi and hi are the vector (transverse) mode, 
and Hij is the tensor (transverse-traceless) mode satisfying 
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b\ t = 0, (3.6) 

ti ti = 0, (3.7) 

H\ = 0, (3.8) 

H ij :j = 0. (3.9) 

Raising and lowering indices of the perturbed quantities are done by 5^ and 5^. 

A. Residual gauge freedom in the Lagrange condition 

The general gauge transformation to the first order is induced by the infinitesimal coor- 
dinate transformation 

x^ = af + ^. (3.10) 
The changes due to the gauge transformation are 

= (3.12) 

In order to leave the Lagrangian condition unchanged, we have 5^ = 0, which leads to 
£",0 = 0. Therefore, 

e = e(x). (3-13) 

For the spatially flat background, £ M are decomposed into each mode, 

e(x) = (T,^L J+ t), (3.14) 

where the vector mode quantity t satisfies = 0. Hence we obtain 

B = B + T(x), (3.15) 

H L = H L --T(x), (3.16) 

(X 

H T = H T -L{x), (3.17) 

h i = h i -£ l (x), (3.18) 
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where an overdot Q denotes d/dt. The vector mode quantity bi and the tensor mode quantity 
Hij are gauge-invariant. 

For a general case of P = P(s), the scale factor takes a complicated form. For simplicity, 
let us take the equation of state as 

P = \e- (3.19) 



B. Scalar perturbations 



The Einstein equations for the scalar perturbations are 



Q-B + H L = 0, 



F t + 4-# L - --V 2 # L + --V 2 H T - --V 2 B = 0, 

a 6 a 6 a 6 a 



H T + 3-H T = \(H L + -{aB) 

a or v «■ 



(3.20) 
(3.21) 
(3.22) 



where V 2 = S^d^dj. We introduce the conformal time rj as 



(it 

drj = — . 

a 



(3.23) 



We expand all quantities in Fourier's series; for instance 



\Er = J d 3 k^ k (r]) exp (ik ■ x) 



(3.24) 



where the subscript k denotes the Fourier coefficient. We define 9 as 

k V 

where k = \k\. We are in position to solve Eqs. (|3.20|) - (|3~^) and obtain 

B k = \ (^-) 2 [*fc(-2 + cos# + 9 sin 9) + $ fc (-2 sinfl + 9 cos 0)1, 
^Lfc = y k 9- 2 (l - cos 0) + $ fc r 2 sin 0, 

^Tfc = (^) 2 (V 1 sin 9 - l - cos - i) + $ fc (V 1 cos + i sin - 9^ 



(3.25) 



(3.26) 
(3.27) 
(3.28) 
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where 6q denotes 9 at the initial time and we used the residual gauge freedom to set initially 
H Tk = 0. 

The initial entropy field s(x, to) is also expressed by the metric. When the initial entropy 
contrast 5 = (s — Sb)/sb is sufficiently small, we obtain from Eq. (|3.1|) 
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Skim) = 3 ^fc(3fl - 2 (cosfl - 1) + 30„ 1 sm6 - - cos# 

+$ fc (-3^ 2 sin^o + 3^ 1 cos^+ 2 sin ^o)]- (3.29) 



C. Vector perturbations 

The Einstein equations for the vector perturbations are 

v2 (^-^ 6 *) =4 0' 6 *' ( 3 - 3 °) 

(cP'hi-ab^' = 0. (3.31) 

Introducing as 

b l (x)=t 1 / 2 V 2 p i (x), (3.32) 

Eq. (|3.30|) is solved to give 

hi = 2(t - to) 1/2 V 2 A(;r) + 4(t - t )" 1/2 A(^), (3.33) 
where we again used the residual gauge freedom (cf. Eq. (|3.18| )) to set hi(x,t ) = 0. 

D. Tensor perturbations 

The equation for the tensor perturbations is 

E i3 + 3-Hij - \v 2 Hij = 0. (3.34) 
a a 2 

This is a homogeneous wave equation in the expanding universe. The solutions are well 



known and we will not discuss the detail here. See, e.g., [11 . 



E. Discussion 

As for the metric, our result is identical to that of the linear perturbation theory (e.g. 
fl3||). However, it should be emphasized our Lagrangian approach does not rely on the 
assumption that the entropy contrast should be small. It is actually an important advantage 
that we can use (or extrapolate) the well-known solutions of the linear theory to express the 
non-linear behavior of the entropy. 

It is worthwhile to mention that there are some choices of a temporal coordinate. Here, 
we have adopted the proper time, for simplicity. We have another choice w M = (h, 0) [|14[| . 
Then, Eq. fl3.4|) is solved as go% = fi(x)/h 2 . However, we find goo = —1/h 2 , which implies 
that it seems rather tedious to solve the Einstein equation. 



IV. CONCLUSION 

The Lagrangian description in the relativistic cosmology has been extended to the case 
of fluid with pressure. It is applicable to even highly nonlinear regions up to the caustic 
formation, since the entropy and the vorticity are obtained exactly along the fluid flow. In 
this approach, dynamical variables are the metric but not material variables thanks to the 
Lagrange condition. The present description includes the dust cosmology as a special case 
(P = 0) §. 

As future subjects along this Lagrangian approach, it would be interesting to study the 



primordial black hole cosmology (e.g. JT5]]) and the averaging problem in inhomogeneous 
cosmologies (e.g. [|l^ -[T8|). Furthermore, self-interacting cold dark matter (field) models 
have been recently discussed as a way to alleviate inconsistencies between the standard cold 
dark matter scenario and the current observation on galactic and subgalactic scales (< few 
Mpc) []T9] pTf| . It may be also an important application of the present formalism to study 



such a nonlinear scalar dynamics, which is modeled by fluids p^ . 
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